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■ Interplay between Kondo effect and trends to antiferromagnetic and spin glass ordering in perfect 

and disordered bipartite Kondo lattices is considered. Ginzburg-Landau equation is derived from the 
£SJ . microscopic effective action written in three mode representation (Kondo screening, antiferromag- 

netic correlations and spin liquid correlations) . The problem of local constraint is resolved by means 
of Popov-Fedotov representation for localized spin operators. It is shown that the Kondo screening 
enhances the trend to a spin liquid crossover and suppresses antiferromagnetic ordering in perfect 
Kondo lattices and spin glass ordering in doped Kondo lattices. The modified Doniach's diagram is 
constructed, and possibilities of going beyond the mean field approximation are discussed. 
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I. INTRODUCTION 



The Kondo lattice (KL) systems are famous with their unusual electronic and magnetic properties, includ- 
ing giant effective, masses observed in thermodynamic and de Haas- van Alphcn measurements^!, unconventional 
superconductivityrl and fascinating variety of magnetic properties!! The overwhelming majority of metallic KL 
systems demonstrate antiferromagnetic (AFM) correlations, .and one can meet all types of AFM order: localized 
spins in U2Zni7, UCdn, Celn3t3; quadrupole ordering in CeBgu; interplay between localized and itinerant excitations 
in several U- and Ce-based compoundsB; puzzling magnetic order of tiny moments in UPt3, URii2Si2, UNI2AI3 that, 
however, rapidly transform in usual large localized moments with dopingtj; quantum phase transition in CeCu6- x Au x 0; 



fluctuation-type dynamical ordering in U(Pti_ x Pd x )3tj; short-range magnetic correlations without long range order 
but with astonishingly wide temperature interval of critical behavior in CeCug and CeRu2Si2Lj. The list is by no 
means exhaustive. The superconducting state in most cases_also coexists with antiferromagnetism, and, apparently, 
jy-} j Cooper pairing in KL is mediated by magnetic fluctuatian^Hij. The dominant contribution to the low-temperature 
\Q thermodynamic is also given by spin degrees of freedomO'El 

On the other hand, all low-temperature characteristics of KL are scaled by a Kondo temperature Tr. These 
characteristics are Fermi-like, but the esaergy scale of "fermion" spectrum is compressed with a factor Tk/ef relative 
to a conventional electron Fermi liquidti Apparently, the AFM correlations due to Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction I partially suppressed by intra-site Kondo effect should be treated as a background for all 
unusual properties of Kondo lattices, so the main theoretical challenge is to find a scenario of crossover from a high 
temperature regime of weak interaction (scattering) between localized spins and conduction electron Fermi liquid to a 
low-temperature strong coupling regime where the spins lose their localized nature and are confined into unconventional 
c"| ' quantum liquid involving spin degrees of freedom of conduction electrons. 

On the phase diagram of disordered KL more exotic possibilities arise, such as non Fermi liquid regime observed, 
O ■ e.g. near the T = quantum critical point in Yi_ x U x Pd3 (see, e.g., Ref. |l3|). In this family of ternary alloys the 
, spin glass (SG) behavior was .discovered in a U concentration range 0.3 < x < 0.5 with a freezing temperature X* 
growing roonotonicallM-.with aJi-l Among other U- based heavy fermion compounds with SG behavior URh2Ge2t-3, 
U2Rh 3 Si5E3, L^PdSisEZI should be mentioned. The effects of "Kondo disorder" were reported for UCu5_ x Pd x in 
Ref. |l8|. Later on the competition between RKKY and Kondo exchange for disordered Ce alloys was discovered 
experimentally (see Refs. [1^ ^fj). The magnetic phase diagram of CeNii_ x Cu x exhibits change of magnetic ordering 
from AFM to FM at x = 0.8, whereas for 0.2 < x < 0.8 the SG state appears only at high temperatures above the 
FM order. Apparently, the Kondo interaction could be considered as the mechanism leading to reduction of magnetic 
moments because increasing Ni contents effectively reduces the strength of indirect exchange interaction, and then, a 
larger temperature stability range of the SG phase appears (see Refs. [l9|- 120| ). 

The competition between the one-site Kondo type correlations and the indirect inter-site exchange is visualized in 
a Doniach's diagram where all possible phase transition and crossover energies are plotted as a function of a "bare" 
coupling parameter a = J/sf that characterizes exchange interaction between the spin and electron subsystem in 
KLH3. Only Kondo screening and RKKY coupling were competing in original Doniach's diagram. Later on it was 
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noticed that the trend to spin liquid (SL) ordering is the third type of correlation- which modifies essentially the 
magnetic phase diagram of KL in a critical region Tjc ~ I of the Doniach's diagramEjE^I. 

In this paper we present a high temperature mean-field description of transitions from a paramagnetic state to 
correlated spin states in KL, which docs not pplate the local constraint for spin-fermion operators. We use the 
Popov- Fedotov representation of spin operatorsc3 to construct the effective action for KL. In this representation the 
local constraint is satisfied automatically. We consider the mutual influence of various order parameters (Kondo, 
AFM, SL and SG correlation functions) and derive a Ginzburg Landau functional (Section II). On the basis of this 
functional we construct generalized Doniach's diagrams that take into account all interplays. The Doniach's diagram 
for a perfect KL is presented in Section III and the influence of Kondo screening on the SG diagram in disordered KL 
is considered in Section IV. 

All existing theories appeal to mean-field approximations that violate the local gauge invariance both in Kondo 
and SL channelEZI. As a result, fictitious second order phase transitions from a free spin (paramagnetic) state to a 
confined spin (Kondo singlet or resonating valence bond SL) state arise in spite of the fact that neither symmetry is 
violated by these transformations. New approach allows us to get rid of assumption about Kondo- type "condensate" 
within a framework of mean field approach. To eliminate fictitious phase transition to a SL state one should refuse 
from a mean field approach to the SL mode. We offer scenarios of continuous crossover from a paramagnetic state 
of localized spins to SL state, where the interplay between critical AFM fluctuations and Kondo screening clouds in 
KL result in " fermionization" of spin excitations at low temperatures (Section V). In section VI the interrelations 
between the theory and real heavy fermion systems is briefly discussed. 



II. DERIVATION OF EFFECTIVE ACTION 



The Hamiltonian of the KL model is given by 

H = E e * c L c ^ + J E ( s & + \ N i n i) w 

ka j ^ ' 

Here the local electron and spin density operators for conduction electrons in a site j are defined as 

jo a 

where f are the Pauli matrices and Cj a = J^ fc c^ a exp(ikj). The SG freezing is possible if additional quenched 
randomness of inter-site exchange Iji between the localized spins arises. This disorder is described by 

H' = J2lji(SjSi)- (3) 

M 

We start with a perfect Kondo lattice. The spin correlations in KL are characterized by two energy scales, i.e., 
I ~ J 2 /ef, and ~ ef exp(— ef/ J) (the inter-site indirect exchange of the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) type and the Kondo binding energy, respectively). At high enough temperatures the localized spins are 
weakly coupled with the electron Fermi sea having the Fermi energy Ep, so that the magnetic response of a rare-earth 
sublattice of KL is of paramagnetic Curie- Weiss type. With decreasing temperature either crossover to a strong- 
coupling Kondo singlet regime occurs at T ~ A# or the phase transition to an AFM state occurs at T = Tjy ~ zl 
where z is a coordination number in KL. If Tn ~ Ak the interference between two trends result in decreasing of 
both characteristic temperatures or in suppressing one of them. As was noticed in Refs Q|8[ in this case the SL 
correlations with characteristic energies A s ~ I may overcome the AFM correlations, and the spin subsystem of KL 
can condense in a SL state yet in a region of weak Kondo coupling. 

To describe all three modes in a unified way one should derive a free energy functional J~(T) in a region of 
T > (Tn, A x , A s ). First, we integrate out the highest energies ~ ef- Here and below we use the dimensionless coupling 
constants J — > J/ef, I — > I/sf, etc. Since we are still in a weak coupling limit of Kondo-type scattering, we may 
confine ourselves by the standard high temperature renormalization of the one site coupling J — * J(T) = 1/ Lti(T/A#-) 
and the second-order equation of perturbation theory in J for RKKY interaction. As a result, one comes to an effective 
Hamiltonian 

H = E ^c\ a c ka + J]T - lJ2 SjSi + ghJ2i-l) j S* (4) 

ka j jl i 
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Here all energies are measured in sf = 1 units, and an infinitesimal staggered magnetic field is introduced that 
respects the symmetry of magnetic bipartite lattice in AFM case (e f is restored in further calculations wherever it is 
necessary). 

To calculate the spin part of free energy J~ S {T) = —T\t\Z s we represent the partition function Z in terms of path 
integral. The spin subsystem is described by means of Popov-Fedotov trickES 

Z s = Tr e-f m = z w Tr e^^'/W) , (5) 

Here (3 = T -1 , N is the number of unit cells, = J2j Nj , and the spin S =1/2 operators are represented as bilinear 
combinations of fermion operators 

S j = (/jt/rt ^ fllfn)/^ Sf = fy jU Sr = f]j jV (6) 
These operators obey the constraint 

•v; £/j CT /,- i. (7) 

In accordance with Ref. ^6|, the Lagrange term with a fixed imaginary chemical potential —iirT/2 is added to the 
Hamiltonian (|l|) . We use the path integral representation for the partition function, 

Z J DcDcDfDf exp A 
Zo J DcDcDfDf exp A 

Then the Euclidean action for the KL is given by 

rfi 

A = Ao- / drHintir) 
Jo 

Ao=A [c,f}= / drV{^v(r)[a r -e(-*V)+/i]c J v(r)+/ 3 v(r)(a r -*7rT/2)/ iff (T)} (9) 

Following the Popov-Fedotov procedure, the imaginary chemical potential is included in discrete Matsubara fre- 
quencies for semi-fermion operators / J(T . As a result the Matsubara frequencies are determined as oj m — 27rT(m + 1/4) 
for spin semi-fermions and e n — 2nT(n+ 1/2) for conduction electrons. In terms of temperature Green's function the 
Euclidian action has the form 

A = A + Aint = ^ C ka G 1 (k)c ka + ^ fic( u n) D Ov{Un)fjo{u n ) 
ka j<7 

+ 2 51 51 C ja {sx)fj 1 -{ui2)fj, r T-{^l)Cj, rr '{£2)S ei -E 2 ,u l -u 2 
jerer' Cm ,0J n 

f 3 ,a> (w 2 )//, 7 (w 3 )^77'/i7'( w 4)^ 1 (10) 

Here the Green's functions (GF) for bare quasiparticles are 

Go{k,ie n ) = - , D%„{iu m ) = — - (11) 

ie n — Ek + H ioJ m — ogh v 1 1 

(y is the index of magnetic sublattice that defines the direction of staggered magnetic field) . 

The first interaction teem in this equation is responsible for low-energy Kondo correlations, and we will treat 
it in conventional mannerE£l. In the RKKY term two .modes should be considered, i.e. the local mode of AFM 
fluctuationsElffid and the nonlocal spin liquid correlationsE2rti3. We adopt for these modes the Neel type antiferromag- 
netism and resonating valence bond (RVB) type spin liquid state, respectively. In accordance with the general path 
integral approach to KL, we first integrate over fast (electron) degrees of freedom. Then in the sf-exchange contribu- 
tion to the action ( |l0| ) we are left with the auxiliary field </> with statistics complementary to that of semi-fermionsEj. 
The spin correlations in the inter-site RKKY term are treated in terms of vector Bose fields Y (AFM mode) and 
scalar field W (spin liquid RVB mode). As a result, Amt is represented by a following expression 
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A mt = -~TW - Tr— Y q Y_ q - Tr- Wp^Wp^ - Tr/, CT ^G (r)^/ iff . (12) 

J let n i — rt n 



'q J qi-q2 



When making a Fourier transformation for non-local spin liquid correlations (the third term in the r.h.s.) we introduced 
the coordinates R = (Rj + R;)/2 and r = Rj — R; for RVB field, so that P, q are the corresponding momenta. Below 
we assume P = and omit it in notations for the SL mode, Woq = Wq. 

Consequent mean- field approach demands introduction of three " condensates" , i.e. three time independent c- fields 
for Kondo coupling, AFM coupling and SL coupling, respectively, that arise as a consequence of a saddle-point 
approximation for all three modes. For .example, the mean-field description of the interplay between the Kondo 
and RVB couplings was presented inE3nd. The undesirable consequence of this approximation is violation of the 
electromagnetic ?7(1) gauge invariance, when the electrical charge is ascribed to initially neutral spin fermion field / 
(see, e.g. Ref. |l2|). According to a scenario offered inE-llEl, there is no necessity in introducing the mean-field saddle 
point for Kondo coupling because the transition to a correlated spin state occurs at T > Tk- In this case the one 
site Kondo correlations suppress the Neel phase transition (reduce — > Tjv) in favor of spin liquid state with a 
characteristic crossover temperature T* > Tn- So we refrain from using the saddle point approximation for the field 
<j> but still use it for the fields Y and W. 

To compactify the equation for the action A we introduce a spinor representation for semi-fermions 

= (ipT fpl ip+QT ip+Ql) > 
and the following definition of the Fourier transform of inverse semi-fermionic Green's function 

/iu m -W p Y$ Y+ \ 

U m [Wp, Yq - + w . (Id) 

1q 1q lU m — " p+Q U 

V Yq ~Yq iu m -W v+ Q J 



The same function in a lattice representation is presented in Appendix I. This operator arises as a result of Hubbard- 
Stratonovich transformation decoupling the magnetic modes Y and the spin-liquid mode W . Then the effective action 
A s acquires the following form 

A s =TrFD- 1 F + A int . (14) 
Now we integrate over semi-fermionic fields and obtain the effective action for a KL model 

A s = Trln (D m \Y, W) + &G o (r)0,) - ^Tr|0| 2 - Tr^Y q Y_ q - Tr-^W qi W q2 . (15) 

J J q ^qi-q2 

Here the argument |</>| 2 appeared in the Green's function D m as a result of integration of the last term in eq. jl^ ) 
over the semi-fermionic fields. 

In a mean-field approximation for two independent modes (neglecting renormalization due to Kondo scattering) 
cq. ( |l5| ) results in a free energy with two local minima reflecting two possible instabilities of high-temperature 
paramagnetic state relative to Neel and SL states. To describe these instabilities one should pick out the classic part 
of Neel field 

Y = {fiNfl^m^S^B, + Y q (16) 

and use the eikonal approximation for the SL field 

W R<r =/A(r)exp(*0) . (17) 

Here M = (Yq) is the staggered magnetization, Y q are the fluctuations around the mean-field magnetization, Q is 
the AFM vector for a given bipartite lattice, e z is the unit vector along the magnetization axis, A(r) is the modulus 
of RVB field, and 6 = (r ■ A(R)) is the phase of this field. 

As is known, in Heisenberg lattices for dimension d > 1, Tjv is higher than the temperature T s i of crossover to the 
SL state, so that the ordered magnetic phase is the Neel phase. Due to Kondo fluctuations that screen dynamically 
local magnetic correlations and slightly enhance the inter-site semi-fermionic correlations, the balance between two 
modes is shifted towards spin liquid phase in a critical region of Doniach's diagram, Tk ~ I- To show this we include 
in the free energy the corresponding corrections induced by the last term in Eq. (H2h. As was mentioned above we 
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refrain form using the mean-field approach to Kondo field, so that the interplay between the Kondo mode and two 
other modes is taken into account by including the Neel mean filed corrections to the semi-fermionic Green's function. 
Then instead of (O) one has the following equation for D^ 1 : 



U m \J* - Kf j ,„ n . _ Ar n • (.18; 




The next steps, i.e. calculations of fluctuation corrections to the stationary point mean field solutions, can be 
performed by introducing the auxiliary self energies 

M(Y, 6) = D~\Y, W) - D m \M, A) 
K^(co ni ,u) n2 ) = -T5^0 3 -(w ni -n)Go(r,n)^(w„ 2 - fi) (19) 

Then the effective action is approximated by the polynomial expansion 

Tr In (D- 1 (Y, W) + X ) = Tr In D" 1 (TV, A) + Tr £ ^-2 [D m (M, A) [M(Y, 0) + } (20) 

n=l 

(the Fourier transform of diagonal part of the Green's function is calculated in Appendix I) . 

Neglecting all fluctuations, i.e. retaining only the first term in the r.h.s. of Eq. ( |20| ) together with quadratic terms 
for AFM and SL modes ([l5|), one obtains the following expression for the free energy per lattice cell: 

0F{N, A) = _ i n [ 2 cosh(/?z/A72)] + - ^ In [2 cosh(/3/ q A)] (21) 

q 

( Iq = —I). The standard self-consistent mean-field equations for order parameters are obtained by differentiating eq. 
( pl| ) from condition of minima of the free energy. These are 

AA = tanh(^Q (22) 

for the Neel parameter and 

„ / r_ A \ 

(23) 



A = -^Kq)tanh^) 



for the real part of the RVB order parameter. Here i/(q) = Iq/Ia- The latter equation was first derived in Ref. [34|. 

Then making the high temperature expansion of Eq. (pl|), one obtains a Ginzburg-Landau (GL) equation in the 
approximation of two independent modes: 

*W A) = + + ^ Trf + c s ;A 4 (24) 

where tn = 1 — Tn /T and r s ; = 1 — T s i/T. The temperatures of two magnetic instabilities are determined as the 
temperatures of sign inversion in the coefficients in quadratic terms of GL expansion Tn — z\I\/2 and T s i — |/|. The 
forth order GL coefficients cn and c s i are positive and depend only on temperature. Up to this point the theory is 
formulated for arbitrary dimension d. In fact, the dimensionality enters the RKKY coupling parameter (see below) 
and determines the number of nearest neighbors z. We consider zl as a universal parameter in further calculations. 



III. DONIACH'S DIAGRAM REVISITED 



To describe the contribution of Kondo scattering to magnetic part of the Doniach's diagram one should integrate 
A over auxiliary field cj> and thus find the Kondo corrections both to Neel and RVB instability points. One should 
consider two cases: (i) Tjy > T s ; (Kondo corrections screen AFM magnetic moments), and (ii) T s i > Tn (Kondo 
corrections enhance nonlocal RVB correlations). 

(i) Kondo screening of AFM order. In this case one takes A = in the Green's function (|l8|). Then adding the 
last term of eq. (|l2|) to the effective action and integrating over the semi-fermionic fields one gets correction to the 
effective action in a form of polarization operators given by the first diagram in Fig. la. 
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D m (I p A) 





MQJ ^^^^ -MqJ 
G m+n (p + q) 



a) 



b) 



FIG. 1. Diagrams for fluctuation contribution to the effective action responsible for Kondo screening corrections to magnetic 
(a) and spin liquid (b) correlations. 

Here the external wavy lines stand for " semi-bosonic" field <j) describing Kondo correlations (seeil). These "semi- 
bosonic" fields are still bosons from the point of view of the permutation relations, but unlike true bosonic fields they 
do not satisfy symmetric boundary conditions, and cannot condense in a state with zero frequency and momentum. So 
the Popov- Fedotov forpialism gives an adequate description of the fact that there is no broken symmetry corresponding 
to Kondo temperatureEJ. The polarization loop is formed by conduction electron propagator Go (solid line) and local 
semi-fermionic Green's function D m given by Eq. ( |l8|) (see Appendix I for the explicit form of these Green's functions). 
As a result the modified effective action is 



q.n 



an (.AO 



(25) 



The logarithmic renormalization of the coupling constant is already taken into account,in J. Therefore the dimen- 
sionless integral 6H includes only contribution due to nonzero magnetic molecular fieldEa, 



<sn(A0 = 



3h(/W) 



- 1 



o 



fAf 2 



\Te F 



(26) 



(see Appendix II for detailed calculations) . This correlation correction should be incorporated in the equation for the 
free energy, so that 



0J : (Af) = f3F{Af, 0) + Trln 



j 



(27) 



Then differentiation of Eq. (p7|) over the Neel order parameter Af, give the following self-consistent equation in the 
vicinity of renormalized transition point, 



A^ = tanh(«) 



a N cosh 2 (/3/ Q AA/2) 
\n(T/T K ) cosh 2 (/3/ Q A0 



(28) 



instead of Eq 

12 



21). Here the Kondo temperature, Tk, is defined as a temperature where the coefficient in front of 
|0n=o| 2 in Eq. (|25|), i.e. the function J -1 — 5H(J\f), turns into zero. It is seen that the screening corrections near the 
Neel transition point are negative, 6H(ftf — ► 0) = —cin{I3N) 2 < 0, so that Kondo screening effectively increases the 
magnetic free energy, and eventually the logarithmic local field corrections reduce the Neel temperature. The results 
of numerical solution of Eq. ( p8[ ) are shown by a bold line in Fig. 2. Inset (a) illustrates the reduction of Tjv in 
comparison with the bare mean field Neel temperature — zepoi 2 /2, where a = J/ef is dimcnsionless coupling 
constant for the Doniach's diagram. 

(ii) Kondo enhancement of SL transition. Now we assume Af = in Eq. (pl|) and subsequent equations. Following 
the same lines as in preceding subsection, one obtains the modified effective action 



-L - OT(J q A) 



Mq)| 2 



(29) 



G 



instead of fl25h, and the polarization integral with the use of the diagram (b) from Fig.l is given by 



OT(/ q A) = 



1 



cosh/3(7 k A) 



1 + 7 k A tanh(/3/ k A) 



k+q 



(30) 



instead of ( p6| ) (see Appendix II). Here £ p = e p — Ep is the dispersion law for conduction electrons near the Fermi 
surface. Inserting the corresponding corrections to the free energy, 



#F(A) = /3 J 7 (0, A) +Trln 



j - 5n(/ q A) 



(31) 



one obtains the corrected self-consistent equation for A. When deriving this equation, the spinon dispersion can be 
neglected since A — * in a critical point. Then one has 



A = -5>(q) 



tanh 



J(q)A 
T 



a si 



iqA 



Tln{T/T K ] 



(32) 



Here a s ; ~ 1 is a numerical coefficient. It is seen from (^) that unlike the case of local magnetic order, Kondo scattering 
favors transition to a SL state, because this scattering means in fact involvement of itinerant electron spin degrees of 
freedom into spinon dynamics. Mathematically, enhancement arises because <5n(/ q A — > 0) = a s ;(/3/ q A) 2 > 0, so that 
Kondo "anti-screening" effectively decreases the magnetic free energy. The results of numerical solution of Eqs. (|3T 
and d32]) are presented by circles in Fig. 2. 



AFM 



SL 




0.05 0.1 0.15 

Coupling constant J/8 F 



FIG. 2. Doniach diagram reconstructed due to Kondo screening (see text for explanation). 



Here filled circles correspond to the region where the AFM order overcomes the SL phase, and the light circles show 
unphysical "suppressed" AFM solutions obtained beyond the region of validity of the mean- field equation (p8|). Two 
other characteristic temperatures, renormalized Tr and T s i, are shown by dashed and solid line, respectively. The 
effects of suppression of Tjv (thin and thick solid lines for bare and renormalized temperatures) and Tk (thin and 
thick dotted lines) are illustrated by the upper and lower inset, respectively. As is seen from the modified Doniach's 
diagram, the interplay between three modes becomes significant in a critical region where the exchange coupling 
constant is close to the point a c = 0.13 where / = Ak in the conventional Doniach's diagram. If the Kondo screening 
is not taken into account, then T^\a) < T^\a) (thin solid and dotted lines in the lower inset). The Kondo screening 
changes this picture radically, and as a result wide enough interval of the values of parameter a just to the right of the 
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critical value a c arises, where the enhanced transition temperature T s i exceeds both the reduced Neel temperature 
Tjy and the Kondo temperature Tk- The calculations of T s i presented in Fig. 2 are performed for d — 2. The similar 
picture exists for d = 3, although the domain of stable SL state is more narrow (for a given value of zl). This means 
that in this region the stable magnetic phase is, in fact, the spin liquid phase. If one descends from high temperatures 
in a hatched region of Doniach's diagram where Tk ~ Tjy, the Kondo scattering suppresses the AFM correlations, 
but the SL correlations quench Kondo processes at some temperature T s i > Tk- As a result the Kondo- type saddle 
point is not realized in the free energy functional in agreement with the assumption used-,above in our derivation of 
Ginzburg-Landau expansion. The preliminary version of this scenario was presented irO. More refined mean field 
approach described here confirms and enhances this scenario, however the SL liquid phase is still described in the mean 
field approximation. Although the local constraint for spin operators is not violated in Popov-Fedotov formalism, the 
gauge phase is still fixedtH, so the next task is consideration of fluctuation backflows described by the higher order 
terms of Ginzburg-Landau expansion. 



IV. ISING SPIN GLASSES IN DONIACH'S DIAGRAM 



In this section we consider the interplay between Kondo scattering and magnetic correlations in case of random 
RKKY interaction ([}]), where the randomness resultsjn formation of a spin glass phase. We consider disorder induced 
by paramagnetic impurities in KL. As was shown inLJ, elastic scattering results in appearance of a random phase S(r) 
in RKKY indirect exchange parameter, 

J 2 \ cos[2p F r - %(d+l) + 5(r)} 



hi = I(r) cs - J ^ rr)d (33, 

where r = \Ri — Rj\, dis the dimensionality of KL . This form of random exchange predetermines two possible scenarios 
of SG ordering. 

(i) Fluctuations take place around a node.jQf RKKY interaction ([33]). This asymptotic is derived from the general 
equation for RKKY exchange parameteiHl^J 

hj = ^ — JZTi ( J (PF r ) 2 [ J d/2-i{PFr)Y d/2 _ 1 {p F r) + J d /2(PFr)Y d/2 (p F r)] . 

(do is the lattice spacing, J v (x) and Y„(x) are the Bessel functions of 1st and 2nd kind). In this case FM and AFM 
bonds enter the partition function on equal footing, and quenched independent random variables Tj can be described 
by Gaussian distribution P(/y) ~ exp(— ifj A r /(2/ 2 ))ll3. The magnetic ordering effects also can be included in our 
approach by introducing nonzero standard deviation AI ^ into the distribution P{hj) that, in turn, results in 
additional competition between SG and— AFM (or, in some cases, FM) states. Recently, the competition between 
AFM and SG regimes was considered inEJ. 

(ii) RKKY exchange fluctuates around some negative value in AFM domain of exchange parameters. In this case 
there is a competition between SG, SL and AFM phases. The third possibility, i.e. fluctuations in FM domain is 
somewhat trivial because in this case Kondo fluctuations cannot significantly change freezing scenario. 

We start with the case (i). To understand the situation qualitatively we make following simplifying approximations. 
First, we consider only Ising-like exchange in the Hamiltonian (J^): 

H> = -J2liiS?S* (34) 

(ij) 

This is usual approximation in thp-jtheory of spin glasses that allows one to forget about quantum dynamics of spin 
variablesEa. In the original papeiu3 the simplifying assumptions (d = oo, separate electron bath for each localized 
spin) were made. Thus the form of spin-spin correlator was predetermined, and these assumptions allowed the 
authors to obtain exact solution in a framework of dynamical mean field theory approach. We refrain from using 
these approximations. -Second, we confine ourselves with the mean field (replica symmetrical) solution of the Edward- 
Anderson (EA) modelc-1 This means that only a pairwise interaction of nearest neighbors is taken into account. 
The number z of nearest neighbors should be big enough (z~ x <C 1) to verify the mean field approximation. We 
consider the interplay between SG and Kondo-type correlations by means of the replica method. We use the approach 
developed in Ref. 45] for the Sherrington-Kirkpatrick modelcj. Both electron and semi-fermion variables are replicated 
(c — > c a , / — > / a , where a — 1, . . . , n), and the number of replicas is tended to zero, so that the free energy per cell is 
given by the limit T — /3 _1 lim„^ (l — {Z n ) av ) / {nN) . Here the replicated partition function is 
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{Z n )av = I] / dlijPihj) II D{clJtJ exp Mole", / Q ] - J* drH lnt (r) 



(35) 



where ^4o (|l0| ) corresponds to noninteracting fermions. 

Averaging over disorder and integrating out high-energy electronic states with the help of replica-dependent 
Hubbard-Stratonovich trick one comes to the following equation 



with 



JJ J Z?{c^/^r}exp^ + ^Tr[X 2 ]+^drTr|0 Q c a 7 a + a / Q c a - j\<f> a \ 



(36) 



Then following the standard pattern of replica theory for spin glassesEHH one fixes the saddle point in spin space 
related to EA order parameter Q. At this stage the initial problem is mapped onto a set of independent Kondo 
scatterers for low energy conduction electrons in external replica dependent effective magnetic field: 



{Z n ) av =^ V (^-- A z{f3I) 2 N{n? + n{n-l)q 2 )+^\n [] J D{f a , <j> a } £ jT exp(^{/ , 0°, y\ x}) ^ 
where J G f(x) denotes dx / V2ir exp(— x 2 /2)f{x), 

A{r, r, y a , x} = ft [(Dtfr 1 ah{ y \ x)] fa-jYl i^i 2 



(37) 



(38) 



and h(y a ,x) — I^fzqx + Iy/ ' z(q — q)y a is effective local magnetic field, which depends on diagonal, and off-diagonal 
elements of Parisi matrix, q = (Sf(O)Sf(t — > oo)) and q = (Sf(0)S^(t — > oo)) (a ^ b), respectively. The latter one is 
the EA order parameter qEA = <Z- Neglecting all fluctuations and retaining only first two terms in the exponent in 
Eq. (|37|), one comes to EA mean field equation for the free energy, 



/^=^[(i-*~) 2 -(i-<z) 2 ] 



ln(2cosh(/3/a; v ^q)). 



(39) 



(see0). Thenjnaking the high temperature expansion, one obtains the Ginzburg-Landau equation in the vicinity of 
SG transition! 



[3T sg = Z<y ^P q 2 T sg - c sg q 3 + d sg q 4 



(40) 



where r sg = 1 — Tf/T and Tf = ^fzl is a spin glass freezing temperature. 

Like in the previous case of ordered KL we incorporate the static replica dependent magnetic field h in semi- 
fermionic Green's functions. As a result, the modified effective action for Kondo fields arises like in Eqs. J25| ) and 



A[y a , x] = In (2 cosh(/3% a , x)) - ^ 



-j-5n(h(y a ,x)) 



Here similarly to Eq. (Eq 



m(h) 



1 



1+0 



( h2 



2 \caah(/3h) J \Te F t 
Finally, performing Gaussian averaging over 4> fields and taking the limit n — > one obtains the free energy 



In (J 2 cosh(/3%, x))/ [1 - JU(0, h(y, x))}) 



(41) 



(42) 



(43) 
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Corrected equations for q and q are determined from conditions dJ 7 (q,q)/dq = 0, dJ r (q,q)/dq — 0. These are 

G 

2cosh(/3h(y,x))/ [1 - jn(0,%,x))]) . (44) 



(45) 



Under the condition i)<la useful approximate equation for C is obtained^: 

cosh 



i rnnl ~ \\ X i ft , 2/- x\ , « 2 (? - ?(1 + 7^ 2 )) . , [ u f uXy/q 
ln(CC{x,q,q)) = -- In (1 + 7 u z (g - g)) + — — — — + In -• 



2 l + 7u 2 (g-g) 



1 + ju 2 (q — q) 



Here the following short-hand notations are used: u — (3Iy/z, C — J/e^ hi (T/T^) and 7 = 2c/ln(T/Tj<-) with 
c = (7r/4 + 2/7r 2 ) ~ 1. We note again that when J = 0, which corresponds to the absence of Kondo interaction, 
C(x,q,q) — 2 exp (^z((3I) 2 (q — q)j cosh(/3 1 x^/zq), and standard EA equation takes place, providing, e.g. an exact 
identity q = 1 

In the vicinity of a freezing point Eq. (Ej) acquires the form 



* of - 

\n(T/T K ) \ln 2 (T/T K ) 



g"=i- ,„^^ r -Q( ; 27^ r ). ( 46 ) 



',1 + 2cz(/3I) 2 ( ? ~- q)/HT/T K ) ) + ° \\n 2 {T/T K ) 



!3IxJ~zq 



As a result of numerical solution of Eqs. ( (46[ ) we obtain the analog of Doniach's diagram for a disordered KL 
where the spi 
respectively) . 



where the spin glass freezing temperatures without and with Kondo screening contributions are shown (T^ and Tf, 




FIG. 3. Doniach's diagram for spin glass transition in disordered KL (see text for explanation). 

Here t) 0) is obtained from GL equation (^) neglecting Kondo screening effect, and Tj was defined from Eqs. difjj ) 

under additional condition d 2 T sg /dq 2 = 0. The influence of Kondo screening on the diagonal element of Parisi matrix 
q is illustrated by the inset (bare value of q = 1 is shown by the dashed line). Like in the case of perfect KL, the 

screening effect is noticeable when Tj, ' ~ . 
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The influence of SG transition on a Kondo temperature for a KL with SG freezing was studied recently in Ref. |49|. 
Although the Kondo effect in this paper is considered in a mean field approximation (i.e. Kondo screening is treated 
as a true phase transition) and a static ansatz was applied for SG, the authors obtained strong reduction of Kondo 
temperature in the same region 



V. CORRELATIONS IN KONDO LATTICE BEYOND THE MEAN FIELD 



The mean-field Doniach's diagram even in its improved form oversimplifies enormously the real picture of interplay 
between three competing modes in effective action ( |l2| ) . First of all, the proximity of three characteristic temperatures, 
Tk, T s i and Tjy means that even when one of them is dominant, i.e. determines the local minimum of the free energy, 
two other define the fluctuations around the saddle point. Second, it is clear physically that only Neel temperature, 
T/v, is a temperature of a real phase transition, whereas T s i and Tk are merely characteristic crossover temperatures. 
The main shortcoming of the mean field approximation is that this approach treats all three modes on equal footing. 
The method described in preceding section allows one to get rid of artificial phase transition at T — Tk, however 
the problems with description of SL phase titiiLyMft., Meanwhile, it is known that the mean-field approximation for 
SL state violates the local gauge invarianc£L^l23'cZJE3 and fixes the phase 9 of SL mode W (|l7|). Second-order phase 
transition from paramagnetic to SL stateEj is an undesirable corollary of this crude approximation, and fluctuation 
corrections to the mean-field solution cannot improve this defect of the theory. 

In this section we consider several scenarios of mode-mode correlations in a system described by the general equation 
(p|),(p^) for effective action A. First, we offer the description of crossover to a SL state, which allows one to bypass 
the mean field saddle point ((2^). It will be demonstrated that the interplay between fluctuations of the fields </> and 
Yq can trigger the transformation of localized critical relaxation AFM modes into SL type correlations without loss 
of criticality. The main idea of our scenarios is that the heavy fermion state of KL is, in fact, unconventional AFM 
state with spin excitations changing their character from Bose-like spin fluctuations or spin waves to Fermi-like spinon 
modes. Next, we consider the behavior of Kondo mode below T s i and describe the quenching of Kondo scattering by 
SL fluctuations in a hatched part of Doniach's diagram (Fig. 2) where the static molecular field is absent. 




Y ' 

j VW\A 



fWW l 




a) 

FIG. 4. Diagrams for fourth and sixth order polarization operators II4 (a) and 116 (b) in effective action responsible for 
mode-mode coupling. 



We demonstrated above that the Kondo screening enhances SL correlations on a level of the mean field approxima- 
tion. Similar effect should exist on a more refined level of interacting fluctuation modes. To find the corresponding 
mechanism we refuse from bilocal representation of spin mode. Instead of introducing the mode W associated with 
gauge non- invariant U(l) field described by the phase 9 in Eq. (|i~7|), we consider the effect of interference of Kondo 
screening modes associated with spins located on different sites of KL. In fact we consider the high-temperature 
precursors of orthogonality catastrophe mentioned by Nozieres in his formulation of "exhaustion problem"EH. In a 
revised scheme we start with the action determined by the Hamiltonian (|l|). Starting with integration over "fast" 
electronic variables (with energies ~ £f), we obtain Ai n t in a form 



Ant = --Tr|0| 2 - Tr-^Y q Y_ q - Tr/ JCT ^G (r)^/ /o 



J 



/II 



4<Pl<Pj 



'q 
TrY, 



zii 



■6<Pl<Pj I! 



(47) 
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Here instead of introducing the scalar mode W we retained higher order terms in Kondo screening fields. These terms 
are illustrated by the diagrams of Fig. 4 

The diagram of Fig. 4a describes interference of Kondo clouds around the sites R, and Zig-zag lines stand 
for AFM vector mode. Like all screening diagrams in Fermi systems it contains Friedel-like oscillating factor. To 
estimate the polarization operator we use the asymptotic form of the electron Green's function in d-dimcnsion at 
large distancesE3H 



G(r,fi) 



ttt CX P 



2e F 



Pft + i Pft 



1 



sgnfi 



Inserting this function in a four-tail diagram of Fig. 4a, one comes to 

1 cos(2p F r- (d+l)f) 



IT 



(2p F r) 



d-l 





' T ' 








) 


\ £ F 


. £ F_ 





(48) 



(49) 



Therefore we expect that this interference correlates with RKKY-type magnetic order, and the interaction between 
the corresponding modes represented by the diagram (b) in Fig. 4 influences the magnetic response in a "critical" 
region of the Doniach's diagram. This response is determined by the fluctuation corrections to Neel effective action, 

5A eff = ] Yl Y a (q,u n )[r 1 ( q )+xoS n ,o]Y a (q,uj n ) (50) 



Here a are Cartesian coordinates, xo = /?/4 is a static Curie susceptibility of isolated spin 1/2 (Fig. 5a). The term in 
square brackets is, in fact, inverse Ornstein-Zernicke correlator ~ &o(Q — °l) 2 + t n a t T > Tn and Q q — > 0. First 
non-vanishing correction to xo is given by Fig. 5b. In this diagram the spins Sj and Sj are screened independently, 
(the wavy lines represent all parquet vertex insertions). In the mean field approach the similar effects are described 
by Eq. @. 




a) b) 

FIG. 5. Diagrams describing local (Curie-type) magnetic susceptibility \o ( a ) an d nonlocal correction taking into account 
Kondo screening of vertices (b). 

Indeed, each vertex correction Ti=jj(u>,e) ~ {4>{e)<f>{e)) gives the jeoptribution ~ l/ln(e/Tx), and integration over 
internal frequency e results in 1/ In (T^/Tr-) correction in Eq. (p8|)c3E3. 

The effects essentially beyond the mean field are described by the diagrams that cannot be cut along a pair of 
electron propagators (solid lines) (see Fig. 6. and 7a). 




a) b) 
FIG. 6. Leading diagrams describing interference of Kondo clouds in magnetic susceptibility (see text for details). 
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The first of these diagrams (Fig. 6a) can be treated as a nonlocal correction to one site spin susceptibility (Fig. 5a) 
induced by interfering flow and counterflow of two Kondo clouds. As a result, the spin-fermion propagator becomes 
nonlocal without introducing the mean field order parameter (|l7|). The next diagram (Fig. 6b) is a kind of "exchange" 
by these clouds in the course of two-spinon propagation. Up to now we exploited the "proximity" effects T > Tk • A 
critical AFM mode given by the Fourier transform of the diagram of Fig. 5a with the wave vector q ~ Q also exists 
in this temperature interval, and, moreover, this mode is dominant in spin susceptibility at T > Tjy. This means that 
the nonlocal contributions of Fig. 6 should be taken also at these q. Due to nonlocality, the temperature dependence 
of spin polarization loop will be slower than the Curie law X/T, and the inverse static susceptibility given by these 
diagrams is 

X q 1 (T) = Xv 1 (T) + x; 1 1 (T) + Iq (51) 

This deviation from Curie law results in delay of Neel phase transition or, in other words, in extension of critical 
regime to temperatures well below TjJ- in accordance with scenario described in Ref. [2^. Magnetic instabilities that 
can emerge at T <C T^ will be the instabilities of spin liquid phase. These instabilities have much in common with 
itinerant fluctuational magnetism considered, e.g. in Refs. p2| , ^3[ 




a) b) 

FIG. 7. (a) Next to parquet approximation for Kondo correction to the magnetic susceptibility; (b) magnetic fluctuation 
correction to single site Kondo scattering. 

The diagram of Fig. 7a with bare spinon propagators gives only local correction to the susceptibility, however 
at T <C Tff where the spinon lines are dressed by the self energies shown in Fig. 6a, this diagrams also becomes 
nonlocal and, therefore contributes in nonlocal term in the r.h.s. of Eq. (fill). The processes taken into account in the 
diagram of Fig. 7b describe the feedback influence of spin fluctuations on the Kondo screening. This diagram together 
with higher order terms of the same type results in dynamical suppression of Tk as a result of appearance of spin 
fluctuation energy uj s f ~ £~ z in the Kondo logarithm. ln(ej?/max{T, u> s t}). This mechanism is effective not too close 
to real TV where the magnetic correlation length £ determining the short-range magnetic order is still comparable 
with the lattice spacing (here z is dynamical critical exponent). 

This schematic description is only scenario of the theory of critical phenomena in KL. We leave discussion of 
fluctuations around SG transitions beyond the scope of this paper. Some details of a new modulated replica symmetry 
breaking schemes, which combine tree and wave-like structures in AFM SG may be found in Ref. |4l]. More detailed 
calculation of critical magnetic and spin glass fluctuations in spin liquid will be published separately. 

VI. CONCLUDING REMARKS 

We derived in this paper the phase diagram for the Kondo lattice model, starting with a high temperature expansion 
of effective action. As a first step, we succeeded in getting rid of one of fictitious saddle points, i.e. we avoided the 
introduction of " Kondo-condensate" averages {ct a fia) used in previous revisions of the Doniach's diagranc3l23. In 
our modified Doniach's diagram (Fig. 3) the renormalized Tk is the lowest of all characteristic temperatures for all 
reasonable values of coupling constant a where one can neglect valence fluctuations. In fact, the mean- field calculations 
of Ref. give similar picture. The feedback of this result is that the strong Kondo regime is unachievable in a critical 
region of Doniach's diagram, and the real role of Kondo screening for small a where Tjv > T s i > Tk is to reduce 
localized magnetic moments and enhance the electronic density of states around Ef- Thus the moderately heavy 
fermion syrStems with relatively big magnetic moments ordered antiferromagnetically arise (Celn3, CeA^) are possible 
examples)Ej. 



13 



In a critical region of Doniach's diagram Kondo screening changes radically the behavior of KL. According to our 
mean field results the conventional AFM order is suppressed at T ~ T s i > Tk- The SL phase that arises instead is, 
nevertheless, close to magnetic instability, and one can expect that spin subsystem eventually orders magnetically. 
If new transition temperature, Tjv is finite, the singlet spinon coupling is incomplete, so that RVBs have residual 
magnetic moments, and these moments are ordered at T — Tk (we emphasize once more that TV marked by light 
circles in a hatched region of the phase diagram of Fig. 3 is not a real transition temperature. It rather designates the 
temperature region where critical AFM fluctuations arise). Of course, the magnitude of these moments is extremely 
small, and one can qualify this type of magnetic order is intermediate between localized and itinerant AFM. In the 
temperature interval TV < T < Tjv the critical AFM relaxation mode characterizes the magnetic response of the 
system. When Tk — 0, one deals with quantum phase transition, and the case Tk < 0, apparently, corresponds to 
short-range correlations existing in a wide temperature interval < T < Tk- This picture describes in gross features 
the magnetic properties of magnetic KL, but any kind of quantitative description will be possible only after realization 
of scenarios for the critical behavior of spin liquid briefly sketched in Section V. 

Now we turn to discussion of conclusions that could be derived from our theory concerning the nature of heavy 
fermion state. The most important one is that the separation of charge and spin degrees of freedom existing in 
KL at high temperatures take place also in a strong coupling regime at T <§; Tk- Indeed, at hight T exceeding all 
characteristic temperatures in KL the spin excitation spectrum is simple structureless peak of the width T around 
zero energy. This peak is manifested as Curie-type magnetic susceptibility and trivial high-temperature corrections 
~ 1/T" to all thermodynamic quantities due to weak paramagnetic spin scattering of conduction electrons, whose 
Fermi liquid continuum exists as independent charge branch of elementary excitations. Since all transformations of 
spin subsystem occur at T > Tk ( at least in a region of a < 0.2 where the valence fluctuations are still negligible), 
this central peak still exist in a strong coupling regime. Below T s i ~ Tr this peak is formed by spin liquid excitations. 
The character of these excitations reminds relaxation modes in a picture of fluctuation itinerant magnetisrr£3'E3 in 
a wide temperature interval down to T co h where the coherent spin liquid regime of Fermi type is established. The 
interaction between SL mode and conduction electrons is the same exchange-type scattering as at high temperatures. 



This coupling constant J is, however, enhanced by the Kondo effect (see Eq. 25). The electrons in a layer of the 
width Tr a*ound Fermi level interact non-adiabatically with spin fermions at low T. As a result the giant Migdal 
effect arisesEI which results in strong electron mass enhancement. So, the heavy fermion state in accordance with this 
picture is a two-component Fermi liquid where the characteristic energies of charge subsystem (slow electrons with 
e < Tk) and spin subsystem (spinons with lj ~ T s {) are nearly the same). 

Exponentially narrow low-energy peak of predominantly spin origin appears practically in all theories of strongly 
correlated electron systems. In the archetypal Hubbard model this peak arises on a dielectric side of Mott-Hubbard 
transition, and still exists on metallic side, where the charge and spin degrees of freedom are already coupled. This 
is the point where the links between Hubbard and Anderson models arise at least on a level of dynamical mean field 
theory (DMFT) valid at d — > ocE3. On the other hand the mean field solution that results in merged charge and 
spin degrees of freedom in a central peak becomes exact in the large- N theories for N = oo saddle pointcZL Recent 
achievements in this direction are connected with confirmation of Noziere's prediction of second scale in Kondo 
latticecil in the limit of exhaustion regime of small electron concentration. At this temperature the "bachelor" spins 
form a coherent Fermi liquid and lose their localized nature. This anticipation was confirmed by recent calculation 
within mean field slave boson approximation of N — > oo theoryE3. In our approach the regime of bachelor spins does 
not arise, because the Kondo coupling remains weak even at T -C Tk (see above), but the spin degrees of freedom 
become coherent at T ~ T co h, so that two coherence scales is the intrinsic property of the model. 

Another aspect of large N theories is the possibility of supersymmetric description that allows combined description 
of spin degrees of freedom in a mixed fermion-boson SU{N) representation^. This approach allowed the authors to 
retain inter-site RKKY interaction in the limit of N — > oo in spite of 1 /N 2 effect of suppression of all inter-site magnetic 
correlations in a standard large N approach. The use of Popov-Fedotov representation allows treatment of different 
magnetic modes described by these operators as " semi- fermions" or "semi-bosons" in different physical situationsES. 
In this paper we appealed to SU (2) symmetry. The general recipe of generation of modes with intermediate statistics 
between Fermi and Bose limiting cases for SU(N) algebra is offered ir£3. In fact the eventual transformation of 
the states with intermediate statistics into true fermions (bosons) occurs only at T — > 0. so this approach may be 
extremely useful for adequate description of quantum phase transitional. 

In principle other collective modes can modify the scenario of AFM phase transition in KL. In particular, the 
low-lying crystal field excitations may intervene the magnetic phase transition in the same fashion as Kondp-clouds 
in our theory. Probably the CeNiSn family of semimetallic Kondo lattices is an example of such interventior£3. 
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APPENDIX I 



To evaluate the contribution of Kondo mode in the expansion ([20|) for effective action, one needs the Fourier 
transform of the Green's function K& (119) . This is 



/ 0„(k)G°(q)0 n (k) 0„(k)G°(q)0„(k + Q) \ 

n (k)G°(q)<^(k) 0„(k)G°(q)0„(k + Q) 

^(k + Q)G°(q)0„(k) 0„(k + Q)G°(q)0„(k + Q) 

\ 0„(k + Q)G°(q)0„(k) 0„(k + Q)G°(q)0 n (k + Q) J 



(AH) 



The components £) m(T (q) of the semi-fermionic Green's function D in (20) are determined by inverting the matrix 
(18). There are normal and anomalous components, 



and 



dTe^(T r / CT (q,T)/ CT (q,0)) 



dre <w » T (T T / <r (q 3 T)/ <r (q+Q,0)) 



iuj m - W a 



(iWm - Wq) 2 - r 2 



{iLu m - w q ) 2 - r 2 ' 



(AI.2) 



(AI.3) 



respectively. Here Y = MIq/2 and W q = 7 q A. 

To perform calculations in real space, one should know the inverse Green's function (n3[) in coordinate representation: 



I lUJm + Yf 

Yf 

V o 





Wji 


\ 


iuj m - Yf 





Wji 









Wij 




iw m + Yf J 



(AI.4) 



It should be noted that the nonlocal term Wji in (AI.4) responsible for SL correlations transforms into diagonal 
term Wq in momentum representation ( |l5| ) , whereas the local staggered field Yj has non-diagonal matrix elements in 
momentum space corresponding to AFM correlations at q = Q. 



APPENDIX II 

The sum of polarization integrals presented in Fig. 1 is given by the following equation 

n„(Y,W q ) = -T £ An.( P )G? n+ „(p + q) 



(AII.1) 



m,er,p 



Only normal component (AI.2) survives in this equation as a result of spin summation. The Neel loop (Fig. la) after 
performing frequency summation acquires the form 



U(Y, 0) = ^2 \ tanh 



4p_ 

2T 



Z P + Y 



X 



cosh(F/T) 
'Y 

tanh ( — 
T 



fe-r) 2 + A 2 (£ P + y) 2 + A 2 
1 1 



(C P -^) 2 + A 2 (£ P + r) 2 + A 2 
C P - Y Cp + Y 



(£ P - y) 2 + a 2 (e P + r) 2 + A 2 



(AII.2) 
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Here £ p 



ip, A = ttT/2. This integral is an even function of the order parameter, 11(1") = II(— Y). Using the 



inequality Y <C ep, two last terms can be simplified, and introducing the integral over the electron band with constant 
density of states po, one has 



il(y,o) = ± Po f* dejtanhr^ 

, PoY 



-CF 



£p 



Y 



2cosh(F/T) e F 



tanh 



T 



Incorporating po in dimensionless variables, one has in the vicinity of Neel point where Y <C T 



n(r,o) = |i in 



/ y 2 



(AII.3) 



(AII.4) 



The logarithmic term is, in fact, included in the renormalized coupling constant J in Eq. fl2q ) for effective action, 
and the remaining terms give Eq. ( p6|) for 8H. Deeper in magnetic phase where 7>T, the Kondo effect is quenched 
by molecular field, so that 



n = In 



o 



•2 



(AII.5) 



Numerical coefficients a,N,bN arising from approximate estimates of the integrals in Eq.( AII.3) are of the order of 
unity. 

The SL loop (Fig. lb) can be estimated for q = 0. After frequency summation it is presented by the following 
integral 



n(o,A) = ~£. 



£ p tanh 



ip 

2T 



'I A 
L, A tanh ' -B- 



X 



T J 2cosh(/ p A/T) 



e 2 + A 2 



(AII.6) 



This function is also even, 11(A) = II(— A), Extracting from (AII.6) the logarithmic term ln(ej?/2T), one comes to 
Eq. d3(i|) for 6H. In a critical region of Doniach's diagram where A <C T, one has 



A 2 



2T 2 ^t=2 +X 2- 
p ^P 



a s i ~ 1. 



(AII.7) 
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